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^ : Abstract 

' Renormalization group (RG) evolution of the neutrino mass matrix may take the value of the 

mixing angle ^13 very close to zero, or make it vanish. On the other hand, starting from ^13 = 
at the high scale it may be possible to generate a non-zero ^13 radiatively. In the most general 
i-^ ' scenario with non- vanishing CP violating Dirac and Majorana phases, we explore the evolution 

pi • in the vicinity of ^13 = 0, in terms of its structure in the complex lA^d, plane. This allows us 

, to explain the apparent singularity in the evolution of the Dirac CP phase 5 at ^13 = 0. We 

also introduce a formalism for calculating the RG evolution of neutrino parameters that uses the 

> 

(<~^ ■ Jarlskog invariant and naturally avoids this singular behaviour. We find that the parameters 

00 ! 

\^ . need to be extremely fine-tuned in order to get exactly vanishing ^13 during evolution. For the 

in 
d 



class of neutrino mass models with ^13 = at the high scale, we calculate the extent to which RG 



QQ ' evolution can generate a nonzero ^13, when the low energy effective theory is the standard model or 
O ■ 

. its minimal supersymmetric extension. We find correlated constraints on ^13, the lightest neutrino 

^ , mass mo, the effective Majorana mass mge measured in the neutrinoless double beta decay, and 

' the supersymmetric parameter tan (5. 
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I. INTRODUCTION 



In the last decade, neutrino experiments have reached a stage where the basic structure of 
the neutrino masses and mixing is more or less clear. We know that the three neutrino flavors 
(z/ct, a G {e,/i, r}) mix to form three neutrino mass eigenstates (i/j,-? G {1,2,3}), which are 



separated by Am, 



m| where, rriij 



denote mass eigenvalues with i,j G {1,2,3}. 



The two sets of eigenstates are connected through Ua = {Up^^ai 
Pontecorvo-Maki-Nakagawa-Sakata neutrino mixing matrix [1 



Ui, where t/pMNS is the 
in the basis where the 



charged lepton mass matrix is assumed to be diagonal. This matrix is parametrized as 



where U is the matrix 
/ 

U = 



PMNS 



e'^^ 
e^^2 



e 



«X3 












\ 



1 



(1) 



C12C13 



S12C13 



S12C23 "~ Ci2S23Si3e C12C23 — Si2S23Si3e 



iS 



■S23C13 



v 



S12S23 ~ Ci2C23Si3e*'^ — C12S23 — Si2C23Si3e*'' C23C13 



(2) 



Here Cij and Sij are the cosines and sines respectively of the mixing angle 6ij, S is the Dirac 
CP violating phase, (pi are the Majorana phases, and Xi are the so-called unphysical phases 
that do not play a role in the phenomenology of neutrino mixing, but whose values may be 
predictable within the context of specific models. The current best-fit values and 3a ranges 
for these parameters are summarized in Table [B It is not known whether the neutrino mass 
ordering is normal (mi < m2 < m3) or inverted (m3 < mi < m2). 









Best fit 


3a range 




[10- 


5eV2] 


7.65 


7.05 - 8.34 


|Amii 


[10- 


-3eV2] 


2.40 


2.07 - 2.75 



sin2 6li2 0.304 0.25 - 0.37 

sin2 6'23 0.50 0.36 - 0.67 

sin2 6'i3 0.01 < 0.056 



TABLE I: The present best-fit values and 3a ranges of oscillation parameters 
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An intriguing situation with the neutrino mixing is that two of the mixing angles, 612 
and 623, are definitely large, while the third angle ^13 is small and may even be zero. Such a 
situation is indicative of some kind of symmetry principle at work. Indeed, there is a whole 
class of models with ^13 ^ that are consistent with data j4|. O23 = 7^/4: and 61^ = are 
allowed by the current data and their origin has been traced to an exact u — r exchange 

n 

symmetry in the neutrino mass matrix [&]. Such symmetries can be realized by models based 



on the discrete non-abelian symmetry groups like A4 [6], D4 



0,^3 B, 



5-4 



9|. Special cases 



of an exact /i — r symmetric matrix corresponding to a symmetry for normal ordering 10| , 
Le — L^ — Lr symmetry for inverted ordering HI] and — Lr symmetry for quasidegenerate 
neutrinos can give ^13 = [12]. Any deviation from this value would indicate breaking of 
these symmetries. Models with discrete abelian symmetries can also make 613 vanish 13]. 
Models involving certain texture zeroes in the neutrino Yukawa matrix or certain scaling 
relations between Majorana matrix elements can also predict zero or almost vanishing ^13 



or those with a 



1( 



SO(IO) models with certain structures for Dirac mass matrices 
SO (3) symmetry can predict ^13 < 10^^ with a normal mass ordering 

Most of the symmetries in these models are obeyed at the high scale, and are broken at the 
low scale by, for example, radiative corrections. If the radiative corrections are large enough, 
any trace of the original symmetry may be wiped out. However in the context of a specific 
model, the compatibility between the high scale symmetry and low scale measurements can 
still be verified. This needs a careful study of the renormalization group (RG) evolution of 
the neutrino mass matrix and the mixing parameters. The basic formalism for calculating 



this evolution has been established in 



17 



like the stability of mixing angles and masses 



points 



33, 



34 



[371, 



25 
351 



38 



scale 



43 



3a 



26 



44 



isl . [lol . I20I ]. Specific features of the evolution, 
21, 13,13, 124], possible occurrence of fixed 



271 ]. evolution of nearly degenerate Majorana neutrinos 



28, 



29|, 



3G, 



31, 



41 



ihe generation of large mixing angles from small angles at the high scale 



41 



421 ] ■ or radiative generation of Ues starting from zero value at high 



45| . have been explored. Threshold effects on masses and mixings, due to 



the decouplin g o 



estimated 



46 



47 



leavy particles involved in the neutrino mass generation, have also been 



48| . These effects can revive 49|, |50| the bimaximal mixing scenario 5l| , 



which predicts 6*13 = 0. 

Analytical expressions for the RG evolution of these parameters have been obtained 
through an expansion in the small parameter ^13 [loj]. For a quantity X G {rrii, 6ij, (pi}, the 



3 



evolution may be written as 

X = Ax + 0(^13) , (3) 

where dot represents the derivative with respect to t = ln(yu/GeV)/(167r^), with /i the 
relevant energy scale. Here Ax is independent of 6'i3, but is a function of 1111,912,023, (pi, S 
in general. In the context of quark-lepton complementarity, approximate but transparent 



analytical expressions were obtained in 52| where a further expansion in the small parameter 
oc + tan^/?) was employed. Here i/r is the Yukawa coupling of the tau lepton and 
tan (3 the ratio of vacuum expectation values of the two Higgses in minimal supersymmetric 
standard model (MSSM). Such an expansion was used to constrain the allowed values of 
mixing angles in the context of tri-bimaximal mixing 5^ and to distinguish between various 



symmetry-based relations at the high scale by comparing the low scale ^13 values 



54| 



A subtle but important^ issue arises in the evolution of the Dirac phase 5 at ^13 = 0. With 
the parametrization in [1^], the evolution formally takes the form 

S = ^ + As + 0(^13) , (4) 

t7l3 

such that the derivative of 6 formally diverges at vanishing ^13, indicating an apparent 
singularity. This is an unphysical singularity: all the elements of the mixing matrix t/pMNS 
evolve continuously, and the peculiar evolution of 6 is related to the fact that 6 is undefined 
at ^13 = 0. This argument is in fact used in [l^ to assert that Ds identically vanishes when 
6'i3 = 0, which leads to a specific value of cot 5 which is a function of {mj,(/)j} at ^13 = 0. 
Ref. [27j has examined this prescription in various limits in the parameter space. 

While the above prescription for choosing the value of S at ^13 = works practically 
when one needs to start with vanishing 613, a few conceptual problems remain. Firstly, when 
6'i3 = 0, the value of 6 chosen should not make a difference to the RG evolution since S is an 
unphysical quantity at this point. Secondly, it is not a priori clear whether the prescription 
would work when ^^13 = is reached during the process of RG evolution. Indeed, getting 
the required value of 6 precisely when ^13 = may seem like fine tuning. The prescription 
in 13], though practical, does not tell us the origin of this apparent coincidence. Here we 
analyze this problem in more detail, and find an explanation in terms of the evolution of 
the complex quantity Ues = sin6'i3e~*'' in the parameter plane Re(We3)^Ini(We3). 

We also evolve an alternative formalism where the singularity does not arise at all. This is 
based on the observation that the set of quantities Vj = {nii, 612, 623, 9l^, (pi, Jqp, Jqp}, where 



Jcp = |si2Ci2S23C23Si3C^3 sm(5 is the Jarlskog invariant and J^p = |si2Ci2S23C23Si3C^3 cos(5, 
have the same information as the set Vs = {w^j, 6*12; ^23? ^13? We therefore write the 

evolution equations in terms of the former set and exphcitly show that the complete evolution 
may be studied without any reference to diverging quantities. We confirm numerically that 
the evolutions with both the parametrizations indeed match with each other and with the 
exact numerical one. 

With the conceptual issue clarified, we numerically study the extent to which ^13 may be 
generated through RG running in the class of models with ^^13 = at the high scale, where 
the low energy effective theory is the standard model (SM) or the minimal supersymmetric 
standard model (MSSM). This evolution turns out to be extremely sensitive to the mass of 
the lightest neutrino mo, the neutrino mass ordering and the Majorana phases. Another ex- 
perimentally observable quantity that depends on these parameters is the effective Majorana 
mass rriee which is explored by the neutrinoless double beta decay experiments. Correlated 
constraints can therefore be obtained on 6'i3, mo and mee, the quantities for which only 
upper bounds are available currently but which may be measured in the next generation 
experiments. For the case of MSSM, it will also depend on the value of tan (3. 

The paper is organized as follows. Sec. [TTl deals with the apparent singularity in the 
evolution of 5. Sec. IIIII calculates the RG evolution in terms of the parameter set Vj. Sec. IIVI 
determines the upper bounds on the value of ^^13 generated through the RG evolution in the 
SM and the MSSM. In Sec. |Vl we summarize our results. 

II. APPARENT SINGULARITY IN 5 AT 6*13 = AND RG EVOLUTION IN THE 
COMPLEX Uez PLANE 

Analytic studies of the evolution of neutrino parameters till date have been mostly per- 
formed with the parameter set V5 = {mj, 6^12, 6^23, 6'i3, 0j, 5}. The RG evolution equations 
obtained are all continuous and non-singular, except the equation for the Dirac CP phase 
5, which is given by 

5 = ^ + As + 0{e,s), (5) 



5 



where 



Ds = % sin 2^12 sin 2^23 X 

mi sin (201 — 5) — {1 + ()'m2 sin (202 — 5) + C'^3 sin 6 



(6) 



A, 



mim2 2 

T 2~'^23 

Am 



sin (201 — 202 



''21 



minis 
Amii 

^2^3 

Am^2 



(44sin(25-20i) + 

S^2 COS 

2^23 sin20i) 

(S12C23 sin (2(5 — 202) + C12 cos 26'23 sin 202) 



(7) 



Here C = ^m\^/ /S.m\2 and C is a constant which depends on the underlying effective theory 
in the energy regime considered. Eq. ([5]) clearly suggests that 5 diverges for ^13 ^ 0. This 
problem is overcome by requiring that Z^^ = at ^13 = 0, which gives the following condition 
on 6 at ^13 = ^]: 

mi cos 201 — (1 + C)"^2 cos 202 — C"^3 



cot 5 



(8) 



mi sin 201 — (1 + C,)m2 sin 202 

The above prescription works for the calculation of evolution when one starts with vanishing 
^13. However on the face of it, it seems to imply that the CP phase 5, which does not have 
any physical meaning at the point 6*13 = 0, should attain a particular value depending on 
the masses and Majorana phases, as given in eq. (IE]). Also, the situation when 6'i3 = 
is reached during the course of the RG evolution has not been studied so far, so it is not 
clear if the prescription needs to be introduced by hand in such a case, or whether the RG 
evolution equations stay valid while passing through 6*13 = 0. Getting the required value of 
5 precisely when 6'i3 = would seem to need fine tuning, unless we are able to figure out 
the origin of this apparent coincidence, and show that this value of 5 is a natural limit of 
the RG evolution. 

The problem also propagates to the evolution of ^13, since it depends in turn on 5: 



013 = Ai, + 0(61, 



A 



13 



sin 2^^i2 sin IQ^ 



23 



^3 
Am2^ 



[mi cos (201 — 5) — (1 + Qm2 cos (202 — S) — (m^ cos 6] 



(9) 



(10) 
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FIG. 1: The left panel shows the evolution in the Re(Z//e3)~Im(^e3) parameter plane, whereas The 
right panel shows the corresponding evolution in the ^13^5 plane. The values of the parameters 
chosen at /^o = 10^^ are: tan/3 = 50, mo = 0.0585 eV^, Am^^ = 4.22 x 10^^ eV^, Amfa = 
3.91 X 10-3 eV^, 6^12 = 32.84°, ^23 = 43.71° and 6^13 = 0.014 rad. The Majorana phases are taken 
to be 01 = 58.9° and 02 = 159.15°. The Dirac CP phase is 124.0° for case A (violet, dash-dotted 
line), 128.447° for case B (red, solid line) and 133.0° for case C (green, dashed line). 

The evolution of all the other quantities, viz. 6*12, ^^23, ""^i, 0i is independent of 5 upto 0{Q\^ 
[iqI . so these quantities do not concern us here. 

In order to understand the nature of the apparent singularity in 5, we explore the RG 
evolution of the complex quantity We3 = sin6'i3e-*'^. We start with three representative 
values of b at the energy scale /^o = 10^^ GeV, with the other parameters chosen such that 
^13 ~ 10^3 at /i ~ 10^ GeV. The left panel of Fig. [1] shows the evolution in the complex 
plane. The right panel shows the corresponding evolution in the ^13-^ plane, with 5 = 27r — 5. 
The following observations may be made from the figures: 

(a) Though all the parameter values at the high scale are very close, and though in all cases 
6*13 decreases to a very small value before it starts to increase, ^13 does not vanish during 
the evolution in all the cases. Indeed, the value of b chosen at the high scale, in order to 
make Q\j, vanish during its evolution, needs to be extremely fine-tuned. This is because 



sin2^i3 = [Re(We3)]' + [Im(We3)]' 



(11) 



so that one needs both the real and imaginary components of Uez to vanish simultaneously, 
which needs a coincidence. Note that when both the CP violating phases 6 and 0j vanish 
at the high scale, Im(We3) = automatically throughout the evolution. Then starting from 



a non-zero value at high scale, ^13 can be made to vanish simply by requiring Re(We3) = 
so that no fine tuning is needed. 

(b) With the definition 5 = 2-k — 5 we have = Si^e^^ and thus 6 is the phase of Ues which 
can be read off easily from the Re(We3)-Im(We3) plot. The values of 6 chosen at fiQ = 10^^ 
GeV are such that 6 is in the third quadrant, so Re(We3) < and Im(We3) < at this scale. 
At the end of the evolution, at /i = 10^ GeV, 6 returns to the third quadrant. During its 
evolution, 6 may change its quadrant zero, one or multiple times. The value of ^13 need 
not vanish completely during the RG evolution, as is represented by the scenarios A and C. 
Scenario B is the one where Re(ZYe3) and lm{Ue3) vanish at the same point, and therefore 
6*13 passes through zero during its evolution. 

(c) In scenario A, since Re(We3) stays negative, 6 simply moves from the third quadrant 
to the second, and then returns to the third in a continuous manner. In scenario C on the 
other hand, 6 has to pass through the fourth, first and second quadrant in sequence to finally 
return to the third quadrant. However its evolution is continuous, the apparent jump at the 
lowest ^13 values in the right panel of Fig. [T]is just the identification of 0° and 360°. 

(d) In scenario B, 6 starts in the third quadrant and moves continuously to the fourth 
quadrant. However it propagates to the second quadrant directly through the origin, thus 
bypassing the first quadrant entirely. Its value at the origin can be well-defined through the 
limit 



cot^o^ lim ^ j.^ '^j"'^^"--^ (12) 

Re(W,3),Im{W,3)^0 Im(We3) Re(We3),Im(We3)^0 ^Im(We3) 

where we have used L'Hospital's rule to compute the limit since both the numerator and 
denominator in this ratio tend to zero at the limiting point. 
Since 

Re(We3) = sin ^13 cos 5 , Im(We3) = — sin 6113 sin 5 , (13) 

we have 

~ Ai3 cos 6 - Ds sin 6 , . 

cot Oq = — : — , (14) 

sm + Ds cos 



and using eqs. ([6]) and (iTOl) . one obtains 



~ micos20i - (1 + C)m2cos202 - C"^3 

COtOo = ; 7 -T ; . (15) 

mi sm 201 — (1 + Qm2 sm 202 



Since 6 = 2n ~ 6, this is equivalent to 



micos20i - (1 + C)m2cos202 - C"^3 

cot do = ] 7 TT ] , (16) 

mi sm 201 — (1 + C)m2 sm 202 



which corresponds exactly to the value of cot S in eq. ([8]), which had been prescribed in 19 1. 
We have thus shown that the prescription follows directly from the procedure of taking the 
limit of S as Re(We3) ^^nd lm{Ue3) go to zero simultaneously. 

The net evolution of 613 and 6 as functions of the energy scale has been shown in the 
top panels of Fig. [2J The evolution of 6 clearly has a discontinuity at ^13 = in scenario 
B, where its value changes by vr. Though the origin of this discontinuity has now been 
well understood, it is important to have a clear evolution of parameters that reflect the 
continuous nature of the evolution of elements of the neutrino mixing matrix ?7pmns- This 
can clearly be achieved by using the parameters ReiUes) and Im(We3)- However, we prefer 
to use the Jarlskog invariant 

Jcp = - sin 6*12 cos 6*12 sin 623 cos 623 sin 6*13 cos^ 613 sin 6 , (17) 
2 

which appears in the probability expressions relevant for the neutrino oscillation experi- 
ments, and is therefore more directly measurable than the real and imaginary parts of Ue3- 
Since Jcp has information only about sin 5, we need its partner 

Jqp = - sin 9i2 cos 612 sin 6*23 cos 6*23 sin ^13 cos^ ^13 cos S , (18) 

to keep track of the quadrant in which 6 lies. The evolutions of {Jqp, Jqp) very similar 
to those of (ReiUes), ImiUes)), as can be seen from the bottom panels of Fig. [2l 



III. RG EVOLUTION EQUATIONS IN TERMS OF THE PARAMETER SET Vj 



We now calculate the RG evolution of the Jarlskog invariant Jcp and its partner J^p as 
defined in f[T5]) . and get to a set of evolution equations that are nonsingular everywhere, 
even at ^13 = 0. The RG evolution equation for Jcp and J^p are obtained as 



with 



JcP = Aj + 0{d,3) , 
j'cp = A'j + 0{9,3) , 



Aj — Cy^S^2^^2^ 23^^23 
— '-^i/T''12'-12''23'-23 



m3 



mi sin 201 — (1 + C)^2 sin 202 
mi cos 201 — (1 + C)"^2 cos 202 — C"^3 
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FIG. 2: The evolution of ^13, 5, Re(Z//e3)! I™(^e3)) Jgv and J^p as functions of the energy scale //, 
in the scenarios A (violet, dash-dotted line), B (red, solid line) and C (green, dashed line). 



We also choose to write the RG evolution for instead of ^13, as is traditionally done. 
This quantity turns out to have a nonsingular behaviour at ^^13 = 0. Moreover, since ^13 > 
by convention, the complete information about ^^13 lies within Q\^. Also, the possible "sign 
problem"^ of ^13 is avoided. In terms of the new parameters Jcp and J^p, the RG evolution 



^ Usually the convention used in defining the elements of ?7pmns is to take the angles Qij to lie in the first 
quadrant. Uez can then take both positive or negative values depending on the choice of the CP phase 
6. In the formulation of eq. (10) the sign of A13 can be such that ^13 can assume negative values during 
the course of evolution and in such situations one will have to talk about the evolution of l^ial- Our 
formulation in terms of Q\^^ as shown in eq. (I23p . naturally avoids this problem. 
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logio (p/GeY) logio (/i/GeV) 

FIG. 3: Comparison of the RG evolution of ^13 and Jcp from the analytic expressions in Vs 
basis (green, dashed line) and Vj basis (blue, dotted line) with the exact numeric one (red, solid 
line). The parameters chosen at the high scale fiQ = 10^^ GeV are: tan/3 = 50, mo = 0.05 eV^, 
Amii = 0.00008 eV^, Am^a = 0.0026 eV^, 6*12 = 34.5°, 6*23 = 42.5° and 613 = 0.5°. The phases 
are taken to be 5 = 40°, (pi = 25° and </>2 = 105°. 

equations for Of^ becomes 

Oh = All + 0{9h) , (23) 
All = ^'^yr^^i'^CP [mi sin 20i - (1 + C)m2 sin 202] 

+ Jqp [nil cos 201 — (1 + ()'m2 cos 202 — C^^s] | • (24) 

Thus the evolution equations in basis Vj are all non-singular and continuous at every point. 
In particular, even when 6 shows a discontinuity, Jcp as well as J^p change in a continuous 
manner. 

In Fig. [3l we show the RG evolution of 9is (left panel) and Jcp (right panel), as obtained 
from the analytic expressions in Vs basis as well as in the Vj basis, along with the exact 
numerical solution, for some chosen values of parameters. It shows that the approximate 
running equations agree with each other to an accuracy of 0(613). 

IV. BOUNDS ON (9i3 AT LOW SCALE 

We now consider all the theories that predict ^13 = at the high scale and try to see the 
nature of running of the masses and mixing parameters with the energy scale. For high scale 
we consider fiQ = 10^^ GeV and implement the symmetry 613 = at this scale, which we 

11 



also take to be the mass of the hghtest heavy particle responsible for the seesaw mechanism. 
We choose this value of /io since it is consistent with the current neutrino mass squared 
differences and seesaw mechanism with Dirac mass of the heaviest neutrino around 1-100 



GeV [55]. This scale is also desirable for successful leptogenesis 56|. However, our results are 
only logarithmically sensitive to this choice and hence our conclusions will be robust against 
variations of /io. Also, this would allow us to compare our bounds with those obtained in 



5J] for specific models like tri-bimaximal mixing at the high scale. The values of the other 
parameters at high energy are chosen such that their low scale values are compatible with 
experiments. For the absolute mass scale of neutrinos, we take the cosmological bound of 



^0 ^ 0.5 eV [57| at the laboratory energy. 

We consider the scenarios where the effective theory below fiQ is the SM or the MSSM. 
We then estimate the maximum value that 613 can gain through radiative corrections. This 
can be obtained from 



^13 = 



to 



(25) 



ICIA 777,3 

^ sin 2012 sin 2923-r-. — — x 



2 ^ "|A7772^| 

|mi cos (201 -6)-{l + C)m2 cos (202 - S) - (rris cos S\ + 0( A^^is, A^) , (26) 

where to = ln(/io/GeV)/(167r2), C = -3/2 for SM and C = 1 for MSSM. Note that we can 
use the parameter set Vs here since apart from the starting point, where 6 is unphysical and 
hence is irrelevant completely, the evolution in terms of this set is also continuous everywhere. 
Moreover it is convenient to talk about Dirac and Majorana phases while putting bounds 
on quantities. In eq. (l26!l . A^^ is defined as 



A«" = In ^ (27) 

327r2 \Mw J Kfi J ^ ^ 

in the SM, where g2 is the SU(2)l gauge coupling, whereas 777,- and Mw are the r lepton 
and W boson masses respectively. In the MSSM, 

1 / \ 2 / 

,MSSM _ / 92mr \ , 2 / /^o 



Numerically, one has A^^ ^ -1.4 x 10"^ and A^^^^ ^ -1.4 x 10~^(1 H-tan^ ^), where tan/? 
can take values upto ~ 50, and so one can treat these quantities as small parameters. We 
explicitly indicate the neglected powers of these parameters in eq. (I26l) . 
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In order to get the maximum ^13 value possible, for any value of the lowest neutrino mass 
mo, all the coefficients of the masses rrii in eq. fl26l) should have the same sign (which we 
choose to be positive) and the maximum possible magnitude. This can be achieved with the 
choice 

20i-5o = O, |202-5o| = vr, (29) 

which gives us 

^max ^ \£\^^^^2e^2sm2923-P^[m, + {l + C)m2+\Cm■sCosSo\] (30) 
< l^ll^ sin2^i2 sin 2^23^^^ hi + (1 + 0^2+ ICI^s] • (31) 

The right hand side of eq. fl3T]) corresponds to choosing the phases shown in Table [IT] for 
eq. (j26ll . As seen, these phases depend only on whether the neutrino mass ordering is normal 
or inverted, and not on the low energy effective theory (SM or MSSM). However, the value 
itself will indeed depend on the effective theory considered. Note that in this procedure of 
bounding 6*13, the actual value of 60 did not need to be used, a considerable simplification 
achieved at the expense of a small overestimation. 
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TABLE II: Phase choices in SM and MSSM that give the maximum radiative correction for ^13. 



To estimate that can be generated at the low scale, we take the optimal values 



of the other quantities in their current 3a allowed ranges 



. We are allowed to do this 



since the corrections to due to the evolutions of the other quantities will formally be 



0{Al) |52|. The quantity that may run quite a bit is 612, however the running is extremely 
small in the SM and 612 always increases in the MSSM, so we use the maximum allowed 
value of sin26'i2 in eq. (pTl) for our estimation. The values of mi, m2 and m3 depend on 
Am2i, Am32, mo as well as the chosen mass ordering. The running of masses and the mass 
squared differences are governed by the Yukawa couphngs of up-type quarks and the U(l)y 
and SU(2)i gauge couplings. For SM, these evolutions depend also on the Higgs boson self 
coupling, and Yukawa couplings of down-type quarks and charged leptons. But 613, as given 

13 



in eq. (l3Ti) . will be independent of these quantities to the leading order in A,- and thus 
considering Am^i, Am32 in the current 3a range is expected to give the correct estimate 
to this order. This assumption can be seen to be valid a posteriori from the comparison 
between analytic and numerical results that follow. 

A. 01 s at the low scale in the SM 
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FIG. 4: Scatter points show the low energy ^13 as a function of the lightest neutrino mass mo at 
the low scale, for both normal (left panel) and inverted (right panel) mass ordering. Each point 
represents a different high energy theory with ^13 = 0. The solid (black) line gives the maximum 
attainable ^13 for a given mo, calculated using the analytic bound in eq. ([3T]) . the current 3a limits 
of the masses and mixings, and the phase values as given in Table HH 



We first consider the case when the effective low energy theory below /io is the SM. 
Running of the masses and mixing parameters is considered from /io = 10^^ GeV to the 
current experimental scale (~ Mz). The scatter points in Fig. [Hare obtained by keeping 
^13 = and varying the other two mixing angles randomly in the range to 7r/2, whereas 
the phases are varied between to 27r. The masses at the high scale are varied within 0-1.0 
eV, so that the lightest neurtino mass mo at the low scale varies between and 0.5 eV. Thus 
each point represents a different high energy theory with ^^13 = at the high scale. The 
upper bound can be analytically estimated through eq. (|3T|) . which depends on the neutrino 
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mass ordering through the phase choices made in Table [TTl and the value of ^ is given in 
eq. m- 

From Fig. H] it is seen that the maximum value gained radiatively by ^13 is rather small, 
being < 3 x 10"'^ in the range < itlq < 0.5 eV for both the mass orderings. Hence if future 
experiments measure 613 greater than this limit, all the theories with ^13 = at the high 
scale and SM as the low energy effective theory will be ruled out completely. If the upper 



limit for mo is brought down by KATRIN 



591 ] to niQ < 0.2 eV, even lower 613 values will 



be excluded for this class of theories. Note that for tuq of this order, the effective electron 
neutrino mass measured by KATRIN will essentially be the same as tuq. 



B. 6^13 at the low scale from MSSM 




FIG. 5: Maximum ^13 obtained at the low scale as a function of the lightest neutrino mass mo at 
the low scale for tan /? = 10, 20, 30 and 50 in the normal (left panel) and inverted (right panel) 
mass ordering. The plots show that simultaneous measurement of ^13 and niQ will help in ruling 
out of a class of high energy theories with ^13 = 0. However there is a strong dependence on the 
upper limit of tan 

When MSSM is the low energy effective theory, the evolution of the neutrino parameters 
is proportional to (1 + tan^/5), as is seen from eq. (!28|) . where tan/3 may take values up 
to ~ 50. Thus, considerably larger running of 613 may be expected at large tan/3. The 
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variation of ^13 as a function of mo is shown in Fig. [51 From the figure it can be concluded 
that with the current hmit of mo, the radiative correction to ^^13 = at the high scale can 
be large enough to reach the present upper bound of 6'i3 at laboratory energy. However, 
for a given mo ^0.1 eV, the maximum ^13 these theories can generate is significantly lower 
for the whole tan j3 range. For example, if mo happens to be 0.08 eV, the maximum ^13 for 
tan/? = 50 is 6*13 ~ 0.12, i.e. sin^ 26*13 ~ 0.056. Such a 6*13 regime will be probed by the next 



generation neutrino oscillation experiments like Double CHOOZ 60|, Daya Bay 6l|], T2K 
62 1 . Since the tritium beta decay experiment KATRIN jsQ] plans to probe mo ~ 0.2 eV 



only, it may not be enough to rule out theories with larger tan (3. 

However, the neutrinoless double beta decay (0i//3/3) experiments will measure the effec- 
tive Majorana mass of the electron neutrino 



mee = |c^2C?3^ie^'^'^'^^^ + sl^cl^ruie^'^^''-^^ + s^3m3| 



(32) 



The value of will allow us to estimate the uiq range, albeit with a large uncertainty 
owing to the complete lack of knowledge of the phases 5, 0i and 02 currently. The present 
upper bound on the average neutrino mass is niee < 1.1 eV 



NEMO 



next generation experiments like COBRA j64|, CUORE |65|, EXO|66|, GERDA [67], Super 



63(1 ■ whereas the proposed 



68( 1 ■ MOON [69] plan to probe niee in the range as low as 0.01 eV < mgg < 0.1 eV. 



Therefore, combined measurement of ^13 and mee may enable us to put some bound on the 
theories with large tan (3. 

The expression for niee in (133) "^^^ be expanded in terms of the parameter 5q = Aml^/rriQ, 
which is small in the range rriee > 0.01 eV, and the small parameter 6^13, to get 



rrie 



mo cos 29 



12 



1 







'12 



^13 



ef.Jml + Ami, + 0{6i, 6^91,, ef,) (33) 



2 cos 2^12 

for normal mass ordering, where the phases are chosen as given in Table [III For inverted 
mass ordering. 



me 



COS 2^12(1 - e) 



m; 



+ \Aml,\ +0{6l,6eel„9 



13^ ' 



where 



6^ 







•-12 



^13 



(34) 



(35) 



2 cos2^i2(l + A) ' "^'V' cos2^i2v/rTA, 
The quantity A = |Am32|/mQ is bounded from below, while for inverted mass ordering ^0 is 
a small parameter (~ O{10~^)) in the range mee > 0.02 eV, so that e is small in this range. 
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FIG. 6: Scatter points show the value of ^13 generated at the low scale as a function of rriee, for 
normal (left panel) and inverted (right panel) ordering. Each point represents a different high 
energy theory with ^13 = 0. Different symbols (colors) correspond to different ranges of tan /3, viz. 
squares (red) for 1.0 < tan/3 < 20.0, diamonds (blue) for 20.0 < tan/3 < 30.0 and circles (gray) for 
30.0 < tan/3 < 50.0. The lines show analytic estimates of 6^^^: solid (red) line for tan/3 = 20.0, 
dashed (blue) line for tan /3 = 30.0 and dot-dashed (gray) line for tan (3 = 50.0. 



The analytic expressions in eqs. fl33|) and flMl) are valid for niee > 0.01 eV and rriee > 0.02 
eV respectively. In this domain of validity, we invert the relations (133!) and (IMl) to obtain 
mo in terms of rriee, and then use eq. fl3T|) for an analytic estimation of Off"^. For the rriee 
values outside the range of validity, one has to estimate numerically the minimum allowed 
rriee for a given mo and then use eq. flHTl) to determine 6^^. These estimations are shown in 
Fig. [6] for various tan/3 values. The scattered points are the low scale predictions calculated 
numerically, which show the correlated constraints in the parameter space of 613 and meg. 
It may be noted that the analytic bounds on ^13 obtained here as a function of rriee are 
generous overestimations, mainly due to the error in the estimation of mo for a given niee- 
Note that bounds on Oi^ at the low scale generated by RG evolution have been studied 
earlier in the context of specific neutrino mixing scenarios at the high scale, like the quark- 
lepton complementarity or tri-bimaximal mixing 54|], or correlated generation of Am^i and 



e 



13 



43 



44| . The bounds obtained in this section, which are applicable not only for all the 
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models with ^13 = at the high scale, but to all the models with ^13 = anytime during 
their RG evolution, subsumes the earlier analyses with specific models. 

V. SUMMARY 

If the neutrino mixing angle ^13 is extremely small, it could point towards some flavor 
symmetry in the lepton sector. There is indeed a large class of theories of neutrino mass 
that predict extremely small or even vanishing ^13. However, such predictions are normally 
valid at the high scale where the masses of the heavy particle responsible for neutrino mass 
generation lie. Below this scale, radiative corrections give rise to RG evolution of the neutrino 
mixing parameters, which in principle can wipe out signatures of such symmetries. In this 
paper, we explore the RG evolution of all such theories collectively. 

The RG evolution with the traditional parameter set Vs = {mi,6ij,(f)i,6} involves an 
apparent singularity in the evolution of the Dirac phase S when ^13 = 0. This singularity 
is unphysical, since all the elements of the neutrino mixing matrix C/pmns are continuous 
at — 0, and in fact the value of 5 there should be immaterial. A practical solution to 
this situation has already been proposed, which involves prescribing a specific value of cot 5 
when one starts the RG evolution of a model with 613 = 0. However, the rationale behind 
this prescription has not been studied before in the context of the nature of this apparent 
singularity. This issue also becomes relevant for the class of models under consideration 
here, since if ^13 is very close to zero at the high scale, it may vanish completely during 
its RG evolution, and getting the required value of 5 exactly at that point seems like fine 
tuning. 

We explore the apparent singularity in 6 by analyzing the evolution of the complex 
quantity Wes, which stays continuous throughout the RG evolution. We find that a fine 
tuning is indeed required, but that is to ensure that ^13 exactly vanishes. In general, if the 
CP violating Dirac and Majorana phases take nontrivial values, one does not pass through 
^13 = even when one starts with ^13 very close to zero. One needs rather finely tuned values 
for the starting values of the neutrino mixing parameters, unless one introduces a symmetry 
like CP conservation, which makes the Dirac and Majorana phases vanish everywhere. Since 
the latter assumption is used commonly in literature, one tends to miss the fact that getting 
^13 — during RG evolution is possible only in a small region of the parameter space. 
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However, if the parameters happen to be tuned such that ^13 vanishes exactly, we 
show that the hmiting value of 5 as ^13 ^ is indeed the one given by the prescription 
mentioned above. We thus put the prescription on a solid footing by deriving it from 
first principles. We also propose an alternate parametrization using the parameter set 
Vj = {iTii, 9i2, 023,01^, (pi, Jcp, J'cp}, where all the parameters are well-defined everywhere 
and any seemingly nonsingular behavior is avoided. 

For models with exactly vanishing ^13 at the high scale, we study the generation of 
nonzero ^13 through radiative corrections. We consider two scenarios, one when the low 
energy effective theory is the SM, and the other where it is the MSSM. The radiatively 
generated ^^13 values are correlated with the absolute neutrino mass scale rriQ. This scale will 
be probed by the future experiments on tritium beta decay, and indirectly by the neutrinoless 
double beta decay experiments. If the value of mo is indeed restricted to the value ~ 0.2 eV 
which KATRIN will probe, the maximum value of 9is generated can only be < 3 x 10~^ in 
the SM scenario. With the MSSM, the running can be much higher for large tan /3, such that 
the current bound of 613 < 0.22 may be reached. In this scenario, we correlate the bound 
on ^13 with the effective neutrino Majorana mass niee to be measured in the next generation 
neutrinoless double beta decay experiments. The whole class of models considered in this 
paper can then be ruled out from future measurements of ^13, rriee and tan 
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